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THE  MOVEMENT  OP  A  FREE  SURFACE  OP  HOMOGENEOUS  GROUND  DURING 

AN  UNDERGROUND  EXPLOSION 
N.  I.  Onis‘ko  and  Ye.  I.  Shemyakln 


The  problem  on  the  propagation  of  seismic  waves  In 
a  seml-unbo\inded  elastic  medium  from  an  unstatlonary, 
spherically  symmetric  source  Is  examined.  This  prob¬ 
lem  was  set  up  In  connection  with  the  study  of  wave 
processes  formed  In  the  grovind  during  underground 
explosions.  The  main  attention  Is  devoted  to  the 
construction  of  theoretical  oscillograms  of  the 
displacement  of  points  on  a  free  surface  of  a  medium. 

It  is  assumed  that  the  wave  source  Is  described  by 
a  fxuiction  chosen  on  the  basis  of  experimental  data 
on  camoxiflaged  explosions. 

§1.1.  Let  a  signal  source  be  set  In  a  homogeneous  isotropic 
elastic  medium  having  Lame's  constants  X  and  |x  and  density  p  on  some 
sphere  of  radius  R^  with  center  at  the  point  o  (Plg.  l) .  In  the  prob¬ 
lem  of  wave  propagation  dvjrlng  underground  explosions,  the  signal  may 
be  described  by  the  velocity  of  a  p6U?tlcle  Up  =  ( t)  dvirlng  Its 

movement  along  the  radius  from  the  point  o  or  by  the  normal  stress 

Or  =  OR(t). 


Pig.  1,  For  stating  the  problem. 
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Let  us  examine  the  propagation  of  a  signal  v^(t) .  Let  us  assume 
that  at  t  <  0  the  medium  Is  at  rest,  and  the  signal  Is  Inserted  at 


t  =  0. 


In  a  problem  with  central  symmetry,  for  an  elastic  medium  the 
movement  of  a  particle  along  the  radius  w  may  be  represented  In  the 


form 


6  r©(£)-|  ,  n-it, 

‘)  =  i7r[-7rj'  — TT’  '' — ~ 


(1.1) 


Here  <t(  4)  and  Vp  are  the  potential  and  velocity  of  the  distribu¬ 
tion  of  longitudinal  waves 


^  at 


(1.2) 


where  Uj^(t)  =  UQf(t)  Is  a  present  function,  u^  the  constant  dimension 
of  the  velocity,  f(t)  some  continuous  function,  and  f(0)  =  0. 

The  equation  for  0(€)  follows  from  (1.2) .  Its  solution  at  zero 
starting  data  ♦(O)  =  ♦»(0)  =0  Is  determined  by  the  function 


(1-5) 


The  displacement  w(r,  t)  Is  determined  by  differentiation  with 
respect  to  R,  and  the  velocity  of  the  particle  Uj^(R,  t)  Is  determined 
by  additional  differentiation  with  respect  to 

In  the  future  It  will  be  useful  to  represent  formula  (1.5)  in  the 

form 


» 

fl,-fico  CO 

\  '-(I')- 5 /«)«-" 

<5,-100  ^  0 


(1.4) 
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Here,  in  the  contour  Intergals  a  >  0  and  >  0. 

The  following  may  be  concluded  on  the  basis  of  (1.4) .  The  solu¬ 
tion  of  the  problem  on  the  propagation  of  a  signal  In  elastic  media 
may  be  examined  first  of  all  for  some  elementary  signal  i) ,  and 
then,  using  (1.4),  one  can  move  to  the  actual  change  of  the  signal  In 
time.  This  corresponds  to  the  replacement  of  the  signal  assigned  to 
the  sphere  R  =  by  a  signal  from  a  spatially  concentrated  source. 

In  the  future  the  role  of  <t^(  ^)  will  be  played  by  the  potentials 


a+ioo 


0>^lCO 


(1.5) 


describing  the  change  In  the  velocity  of  the  particles  In  time  for  a 
concentrated  source  of  type  6(t)  and  the  change  of  the  point  source 
In  time  as  e(t)  respectively.  In  these  formulas  6(t)  Is  Dirac's  func¬ 
tion,  and  e(t)  Is  a  step  function 

,l  ((>0)'  -  — 


Several  transformations  are  necessary  In  order  to  use  Formula 
(1.5)  In  the  cylindrical  coordinate  system  (r,0,z)  .  As  a  result,  the 
first  formula  of  (1.5)  takes  the  form 


Q  (»)  '  ^ 


(1.6) 


Here  h  Is  the  ordinate  of  the  center  of  the  source,  £  the  ordin¬ 
ate  of  the  point  of  observation,  and  (1)  the  contour; 


S  =  — n]/r*+ (/j  — z)*,  t' =  /  +  n/fo, 


(Rot  »  G>  0) . 


Formula  (1.6)  describes  the  disturbances  emanating  from  the  sovirce 
In  the  direction  of  the  negative  semlaxls  z. 
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In  conclusion  to  Section  1,  let  us  Introduce  formulas  analogous 
to  (1.4)  for  a  signal  of  the  second  type.  On  a  sphere  of  radius  R*Rq, 
the  normal  stress  Oj^(t)  =  —  assigned,  where  Is  the 

dimensionality  constant  of  the  stress,  g(t)  Is  some  continuous  func¬ 
tion,  and  g(0)  =0. 

In  this  case.  Instead  of  (l.5)>  we  have  the  formula  [l] 


(1.7) 


and  Instead  of  (1.4) 


(1.8) 


Where 

00 

f'W-HT  S 

d-to>  ^ 

The  e35)resslon>  for  the  potentials  of  elementary  signals  In 
spherical  and  cylindrical  coordinates  preserve  the  previous  form. 

Let  us  derive  formulas  for  calculating  the  main  normal  stresses 


'^^(R^  t)  and  a^(R,  t) ,  at  an  arbitrary  point  In  the  medium,  which 
correspond  to  the  source  o_(t) 

A 


0  ®*(5)  ,  i-r  <Pi(|)  .  „._Ji 

- -p - - IF-’  -  p 

0  V\<1»'(5)  o ®JiL 

— 9 - — -77  ^  * 

(1.9) 

Formulas  (I.9)  were  derived  using  (l.l),  (1.7)  and  Hooke's  law 
and  describe  the  complex  change  In  the  stresses  with  distance  as  com¬ 
pared  with  the  change  In  pressure  In  an  analogous  problem  for  a  fluid. 
As  follows  from  (1.9)^  the  presence  of  second  and  third  terms  In  these 
formulas  Is  fully  explained  by  the  resistance  of  the  medium  to  dls- 
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placement  (|i,  0)  .  If  p,  =  0,  l.e.,  *  0  (the  propagation  rate  of 

transverse  waves  Vg  vanishes  at  p,  =  0),  then  from  (I.9)  follows  the 
solution  of  the  problem  for  a  fluid 

(1.10) 

The  first  terms  In  (1.9)  determine  the  asymptotic  behavior  of  the 
stresses  at  R  -»  00  .  If  the  concept  of  the  length  of  a  nonstationary 
wave  X  *  VpT  Is  Introduced,  where  T  Is  some  characteristic  time  of  the 
signal,  then  the  first  terms  In  (I.9)  will  describe  the  behavior  of  the 
stresses  In  time  with  high  accuracy  at  distances  R  »X.  We  are  mainly 
Interested  studying  the  disturbance  fields  at  distances  con^parable  to 
the  wave  length.  In  this  case  we  cannot  be  limited  to  asymptotic  formu¬ 
las  and  should  tvirn  to  accurate  formulas. 

2.  Let  \is  examine  the  problem  for  an  elastic  half -space  z  ^  0, 

If  within  It  at  the  point  z  =  h  an  elementary  signal  Is  sent  which 
describes  a  longltudlnal-wave  field  having  a  potential  of  the  type 
(1.5)  for  a  signal  of  the  type  v>j^(t)  . 

The  solution  can  be  set  up  either  by  the  method  of  V.  I.  Smirnov 
and  S.  L.  Sobolev  [2]  or  by  the  method  of  partial  division  of  the 
variables  [5] . 

As  Is  known,  when  a  wave  descrlbabl*  by  a  longitudinal  (or  trans- ' 
verse)  potential  ♦q(|) strikes  a  free  surface  z  =  0;  waves  descrlbable 
by  longitudinal  <p^{r,z,t)  and  transverse  '('j^(r,z,t)  potentials  are  re¬ 
flected  the  surface.  The  potentials  of  the  reflected  field  may  be  re¬ 
presented  by  the  following  expressions; 

00 

«p,(r,  z,  Jiiiz,  t,  l{)Jo(kr)dk 

n 

CO 

z,  0  =  I,  k)Jt(kr)dk 
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«+ico 


«— ‘ICO 

‘)“4r  \  >',(*. 


Here 

.V,  (*,  5)  =  -  ~  Xo  (k,  s),  Y,{k,  s)  =  -  ^  A'„  (A-, 

ot.  *  /rtV  +  /i="  p,  .  yUW^fW,  g  =  2A*  + 

^  —  4*%tjPi,  r  =.  i'*  + 


(1.11) 


In  this  and  are  represented  by  Melllnovskly  contour  Inte¬ 
grals,  and  X^Ck,  s)  should  be  taken  according  to  (1.6)  In  the  form 


or  according  to  (1.5)  In  the  form 


(i.ia) 


x.(&.  = 

(1.15) 

On  the  basis  of  Functions  (l.ll),  (1.12)  and  1.15),  formulas 
determining  the  reflected  displacement  field  for  a  point  source  of 
the  type  (1.12)  In  the  half-space  z  ^  0  may  be  reduced  to  the  form 


Wrir,  2,  ^ap 

^  ^  ry  —  AaP  _ 


Wt 


y  f  1  ry  _  /.a3 

(r,  2,  0  =  —  l^HT  U6*+ 

-  y^«-*"^'‘“l  ^ 


(1.14) 

Here 

,  =  Kr+TV.  r  =  7-.  :i  =  vTT?.  &  =  2  +  ;* 

p 
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Is  Melllnoyskly's  contour  ReC*a  >  0;  arg  «  »  arg  p  «  0  at 

C  >  0. 

The  displacement  field  described  by  Formulas  (l.l4)  Is  caused  by 
movements  In  a  longitudinal  wave  from  an  elementary  source  (1.12) 

^  ^  [fl6  (/  -  aJ{)  -f  1] 

(1.15) 

The  entire  field  of  displacements  (l.l4)  and  (1.15)>  correspond¬ 
ing  to  source  (1.12),  will  be  studied  In  Section  3  and  Is  subsidiary 
In  this  problem. 

The  main  attention  will  be  devoted  to  the  study  of  the  displace¬ 
ment  field  on  the  surface  of  an  elastic  medium  In  order  to  obtain  ac- 
cvirate  oscillograms  of  the  motion  Of  points  on  the  surface  for  signal 
(1.15).  In  this  Expressions  w_(r,  0,  t)  and  w„(r,  0,  t) ,  determined 
below,  describe  the  total  field  of  movement 


«;„(r.O.O-2A5  A(*r) 

0  (i> 

(r.  0. 0  -  -  ^  W.  (*0  (i  ^  «p  [*  (c  -f 

0  (0 

^A  ■»  2i<«A«') 


(1.16) 


In  conclusion  to  Section  2,  let  us  note  that  the  solutions  of 
(l.l4)  and  (l.l6)  obtained  by  the  method  of  partial  division  of  the 
variables  permits  an  elementary  Investigation  of  the  asymptotic  be¬ 
havior  of  the  components  of  the  displacement  field.  This  was  demon¬ 
strated  by  the  works  of  a  Leningrad  group  under  the  direction  of  O.I. 
Petrashen* . 

In  this  problem,  the  asymptotic  study  of  the  formulas  obtained 
will  be  subsidiary  to  the  qualitative  study  of  the  properties  of  the 
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field  on  the  basis  of  Formulas  (2.4)  and  (2.5)  . 

In  Section  3  we  shall  examine  the  qualitative  wave  pattern  in  an 
elastic  half-space.  Detailed  investigations  of  the  accurate  solutions 
of  problems  on  the  dynamic  theory  of  elasticity,  carried  out  by  a  group 
of  Leningrad  dynamics  specialists  [3,4  and  l],  allow,  when  explaining 
the  qualitative  wave  pattern,  the  fundamental  properties  of  the  wave 
fields  to  be  established  quickly  in  concrete  problems  if  solutions  of 
the  type  ( 1 .  l4)  «ire  set  up . 

These  investigations  were  made  using  asymptotic  methods  for  con¬ 
tour  Integrals  (the  method  of  stationary  phase)  with  subsequent  inte¬ 
gration  according  to  Fourier — Bessel  and  are  Intended,  chiefly,  for 
studying  the  field  at  long  distances  from  the  source  (R  >>  X)  .  In 
this  a  single-valued  correspondence  is  set  up  between  the  properties 
of  the  Integrands  of  the  contour  Integrals  from  (l.l4)  and  the  quali¬ 
tative  properties  of  the  wave  field  [5]. 

3.  Without  repeating  the  ciunbersome  calculations  of  the  princi¬ 
pal  terms  in  the  asymptotic  representations  of  the  disturbances  in  the 
vicinity  of  the  wave  front,  let  us  Introduce  the  picture  of  the  fronts 
of  the  wave  field  at  some  moment  in  time  (Fig.  1.)  emd  write  the  cor¬ 
responding  approximate  formulas. 

1®.  In  the  vicinity  of  the  front  of  a  direct  longitudinal  wave 
Pq,  the  displacement  field  is  described  by  Formulas  (1.15)  for  the 
source  (1.12)  . 

If  Formula  (1.4)  is  taken  into  account  and  we  move  to  a  given 
change  of  the  source  in  time,  the  main  parts  of  the  disturbances  at 
V^r^  +  (h  —  z) ^  »  X  will  have  the  form 


//V  • 


x/Al), 


Wt  ^ - — ,T-  I  i  —  x7,  (;)  (X  =  sin  a) 


(i.n) 
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where  a  Is  the  angle  of  Incidence  of  the  wave  Pq  to  the  boundary. 

If  (1.17)  Is  examined  along  some  plane  z  =  0,  which  passes 
through  an  \inbovinded  elastic  medium  at  a  distance  h  from  the 

center  of  the  soxirce,  then  w^,  decrease  as  in  an  unbounded  medium: 

-1  -P 

R  at  R  CD,  and  w.  will  decrease  as  R 

z 

2®.  In  the  vicinity  of  the  front  of  a  reflected  longitudinal 
wave  the  main  terms  In  the  approximate  asymptotic  formulas 
have  the  form 


(g) 


/.  (6.), 


(j:)  y  I  —  j:- 
»■,.  V tj- 


/l (fel) 


(1.18) 


Here 


G. 


(2tV  -  I )?  -  V 1  - 

—  I )■  +  4'r‘*'*  Y '  —  V 1 
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describes  the  change  In  amplitude  of  the  reflected  wave  as  a  function 
of  the  angle  of  reflection  a.  It  Is  Interesting  to  note  that  the  entry 
sign  of  the  reflected  wave  can  change,  as  the  graph  of  the  function 
Gjl  for  7  =0.6  attests  (Fig.  2)  .  In  accordance  with  these  graphs  In 
Pig.  2,  the  arrow  on  the  front  of  the  wave  Indicates  the  direction 
of  the  motion  of  the  particle  when  entering  this  wave. 

Formulas  (l.l8)  Indicate  the  nattire  of  the  change  of  movement  In 
time,  as  In  the  direct  wave  Pq. 

3®.  In  the  vicinity  of  the  front  of  a  transverse  wave  formed 
together  with  the  reflected  wave  P^,  the  corresponding  approximate 
formulas  have  the  form 


V‘ 

(:,(»)/,  0-4) 


(1.19) 


Here 


<M,  *5  4r^- VT- -  2Tr-x^) 


and  t  j  Is  the  arrival  time  of  the  transverse  wave  at  an  observation 
polnt^  with  coordinates  (r,  z)  such  that  =  ^2  +  ^2* 
the  running  time  of  the  wave  Pq  to  the  boundary  z  =  0,  and  tg  Is  the 
running  time  of  the  transverse  reflected  wave  S^.  Let  us  recall  the 
mechanism  of  the  formation  of  the  transverse  wave  S^.  At  each  point 
on  the  surface  z  =  0  at  which  the  longitudinal  wave  Pq  arrived,  com¬ 
pressions,  and  displacements,  which  are  then  propagated  with  different 
velocities  Vp  and  Vg,  Pig.  1  shows  elementary  waves  Sj  leaving  a  free 
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surface  and  going  Into  the  depths  of  space.  The  transverse  wave 

will  be  the  envelope  of  these  waves^  the  eqviatlon  of  this  envelope 

at  some  fixed  moment  t  =  t  . 

f 

The  function  for  all  x  and  'i  will  remain  positive  (Fig.  3)  t 

this  attests  to  the  fact  that  the  entry  signs  of  the  transverse  wave 
at  the  moment  t  =  t^  In  the  notations  w^(r,  z,  t)  and  w^(r,  z,  t)  at 
all  distances  r  at  fixed  values  of  £  and  h  are  not  changed:  In  entry 
of  a  transverse  wave  a  particle  Is  moved  above  and  ediead  of  the  source, 
moving  along  the  tangent  to  the  front  of  wave  (Plg.  1) . 


Pig.  2.  Graphs  of  G  (x)  for  7=  Pig.  J>.  Graphs  of  G^Cx)  . 

*0.4 — 1,7=  0.5 — II  ^  and7=  0.6 — III.  Values  of  7  same  as  In  Pig.  2. 

As  was  noted  above,  for  the  problem  at  hand  a  correspondence  Is 
set  up  between  the  properties  of  the  Integrands  of  the  contour  Inte¬ 
grals  and  the  properties  of  the  disturbance  field  In  the  solution  of 
(l.l4) .  In  the  Mellln  Integrals  from(l.l4)  the  Integrands  have  simple 
poles  at  the  points  f=tl^(0<<^  <l)  of  the  complex  plane  f.  The 
residues  at  the  points  C  =  1  1  and  subsequent  elementary  Integration 
with  respect  to  k  allows  It  to  be  established  that  the  property  of  the 

^The  time  tf  Is  determined  by  the  formula 

*  *’p  "» 

Here  for  fixed  values  of  r,  z,  h  Is  determined  from  Snell «s  law 
7  sin  a*  sin  p  ,  where  J  Is  the  angle  of  reflection  of  a  transverse 
wave 
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disturbance  field  connected  with  this  Is  none  other  than  a  surface 
wave  of  the  Rayleigh  type  [4].  In  Section  2  this  disturbance  was  ex¬ 
amined  In  more  detail,  by  accurate  formulas  for  the  entire  disturbance 
field.  It  should  be  noted  that  the  principal  part  of  this  disturbance 
Is  propagated  with  a  rate  v_  =  ^  <  v  and,  therefore,  all  pf  this 

disturbance  appears  mainly  during  the  observation  times  t  z  r/v  for 
given  r  >  h  and  r  >2. 

This  observation  was  made  by  analogy  with  known  results,  but  may 
be  easily  demonstrated  by  calculating  the  residues  at  the  points 
^  =  t  1  ♦  In  (l.l4),  and  by  careful  analysis  of  the  properties  of  this 
wave  In  ovu?  problem.  Pig.  1  shows  this  by  a  conventional  symbol  and 
the  letter  R  Is  the  surface  wave. 

Finally,  let  us  Introduce  the  approximate  formulas  for  the  com¬ 
ponents  of  the  displacement  field  on  the  free  surface  of  the  medlxan. 

At  points  on  the  surface  2=0  these  formulas  describe  the  total  field 
of  movements  Pq  +  and  may  be  obtained  either  by  asymptotic 

methods  directly  from  (l.l4)  or  by  summation  of  the  main  parts  of  (1.15)* 

(1.18)  and  (1.19) 


Wr(r,  0,  0~ 

Ut(r,  0,  0»- 


'iHoffn-Cj  (x) 


/id) 


(x) 


P„Vr>  +  /,» 


Ml) 


(1.20) 


Here  {  Is  according  to  (1.20);  the  fvmctlons 

=  rfe  T'  Ki-tV  .  ^^4  =  ^  (1  -  2r-x^) 

A (x)  =  (2t*x»  —  1)*  +  4t*xS |M  —  X*  Vl—yV 

As  follows  from  (l.20),  the  components  of  the  movement  field  de¬ 
crease  differently  with  distance,  as  Is  Illustrated  by  the  graphs  in 
Pigs.  4  and  5. 
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In  connection  with  the  graphs  of  C}j(x)  and  It  should  be 

noted  that  at  r  »  X  the  laws  of  amplitude  reduction  with  distance  for 
the  displacement  field  and  rate  field  are  the  samej  It  vetrles  only  as 
a  curve. 

These  qualitative  properties  of  wave  disturbances  allow  the  funda¬ 
mental  types  of  waves  In  a  medltam,  their  kinetic  and  dynamic  features, 
to  be  shown.  As  usual,  the  wave  propagation  rates  and  arrival  time  of 
waves  at  a  given  point  of  observation  pertain  to  the  first.  The  asymp¬ 
totic  laws  of  the  damping  of  amplitude,  the  comparative  Intensities  of 
the  manifestation  of  various  waves,  and  their  entry  signs  on  a  theoretl 
cal  oscillogram  refer  to  the  second. 


At  the  same  time  It  Is  almost  apparent  that  these  formulas  are 
hardly  applicable  for  quantitative  comparisons  with  experimental  data, 
save  for  a  medium  which  Is  very  nearly  elastic,  eind  at  dlst£uices 
r  »  A  .  In  order  to  obtain  more  accvirate  descriptions  of  elastic  wave 
disturbances  In  the  vicinity  of  the  fronts  of  volume  waves  It  Is  neces¬ 
sary  to  either  calculate  the  successive  terms  of  the  asymptotic  expan- 
stlon  or  to  calculate  by  accurate  formulas. 

§2.1  Accurate  Investigations  of  the  displacement  field  may  be 
carried  out  either  by  solutions  obtained  by  the  Smirnov — Sobolev 
method  [2]  or  by  formulas  of  type  (l.l4)  or  (l.l6) .  Bearing  In  mind 
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the  problem^  let  us  limit  ovirselves  below  to  the  study  of  the  displace¬ 
ment  field  on  the  surface  z  =  0  (l.l6) . 

Ogurtsov  and  Petrashen*  [4]  have  Indicated  a  simple  method  re¬ 
ducing  and  formulas  of  the  type  (l.l6)  to  real  Integrals.  In  this. 

In  the  contour  Integrals,  the  location  of  segments  along  the  Imaginary 
axis  from  the  lower  point  of  branching  to  Infinity  are  chosen  from 
(l.l6)  (symmetrical  to  the  upper  and  lower  points  of  bramchlng  to  In¬ 
finity),  and  the  contour  Is  deformed  to  the  edges  of  the  segments. 

The  subsequent  rearrangement  of  the  order  of  Integration  and  In¬ 
tegration  with  respect  to  k  allow  (l.l6)  to  be  represented  In  the  form 
of  the  following  simple  expressions  In  the  dimensionless  vcu'lables  (  , 
e  ,  t* : 


“'m  =  4  [Ur,  +  Urn  +  Ur>.],  A'  =  ^  , 
~  7*  jH  *+•  u y,,],  *  =5  -JL. , 


Here 


Ur„  =  r-I 


i 


"  W  ^  «-■  i) 

P  =  V\v  +  (1  -  r’O-)  -  -f  (1  -  t={)») 

-  5’ -'!*<' -7=0=)  t)  ’  ®  =  Iy 

rfi  =  4  (2  -  -  T*- 


(2.1) 


and  f  Is  the  root  of  the  equation 


(2  _  0)=  ^  /,  j/rry* 

1 _ ?i 


«  -  7*  (I*  -f.  • 


•  ff  2  —  0*  .  tf) 

‘'-“.TTT'""! 
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The  Items  and  are  represented  by  the  real  integrals 

Uy  ..  r  /  -  m. 


Here 


4  t  (2^y£V^(„yT=pi5sinf +rXo<»|!)St 


« *f'  (2- xnvx*— 1  yi— T*>-*  ■!_  Ti  yA- 

Va-—,-\ - F(S5  *■”  2  yf; 


P(k)  =  (2-  A*)*  -1-  16  (i*  -  1)  (1  -  t“X»), 


(2.2) 


and  Pj^  and  are  obtained  from  pand  «P  by  substitution  of  \  for 

These  formulas  make  it  possible  to  conduct  accxirate  studies  of 
the  displacement  field  on  the  surface  of  an  elastic  medium.  Prom 
physical  considerations  it  is  clear  that  the  calculation  should  be 
made  atj^+n^^  1.  After  calculating  (2.1)  and  (2.2),  transition 
to  the  soxirce,  which  varies  in  time  according  to  the  law  f(t),  is 
accomplished  using  the  formula 

w(r,  a  <)  =  (/,(/  — T)H»,(r,  s,  t)dx 

*  (2.5) 

where  , 

A  (t)  =  /  (t)  —  ^  /  (r)  exp  (r  —  t)  rfv 
0 

Let  us  make  two  comments  in  connection  with  formulas  (2.l),  (2.2) 
and  (2.3) .  The  solution  for  the  displacement  field  in  the  elastic 
problem  can  be  written  as  a  function  of  the  dimensionless  variables 
i  t*.  This  means  that  for  all  elastic  media  having  various  par¬ 
ameters  Vpj  Vg,  P  ,  but  having  fixed  7 ,  in  the  problem  on  may  calculate 
the  \inlversal  curves  of  w^(t*)  and  w^Ct")  at  a  given  point  of  observa¬ 
tion.  In  this  a  scale  in  units  of  k  is  chosen  for  the  ordinate  and 
the  dimensionless  time  t*  is  plotted  along  the  abscissa.  In  this 
problem  the  notation  of  the  velocity  of  particles  on  a  sphere  of 
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radius  Rq  Is  chosenfas  the  signal  f(t) .  If  the  signal  parameters  Uq 
and  Rq  are  flxed^  during  transition  from  one  elastic  medium  to  another 
the  scales  of  the  displacements  will  vary  In  Inverse  proportion  to 
the  velocity  of  the  transverse  waves. 

If  (2.5)  or  (2.1)  Is  differentiated  with  respect  to  time  (moving 
I o  theoretical  oscillograms  for  the  velocity  of  particles  on  the 
svirface  of  the  groimd) ,  the  scale  factor  for  particle  velocity  may 
be  established: 


(2.4) 

If  It  Is  assvimed  that  Uq  and  Rq  are  fixed  values  for  various 

elastic  media,  then  It  can  be  demonstrated  that  at  equal  dimensionless 

distances  the  amplitudes  of  the  particle  velocities  vary  In  Inverse 

proportion  to  1  ,  The  ratio  7  =  v/v  for  most  real  ground  [6]  varies 

s  p 

little:  Y  =  0.4  to  0.6,  therefore,  the  corresponding  amplitudes  of 
the  particle  velocity  will  differ  little. 


Pig.  6 


^  Rq  should  be  taken  as  that  distance  from  the  charge  at  which  the 
effect  of  the  free  surface  upon  the  recording  of  f(t)  Is  not  substan¬ 
tial. 
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These  remarks  also  determine  the  direction  of  the  numerical  in¬ 
vestigations  whose  res\ilts  are  given  below. 

2.  Let  us  tedce  as  a  basis  a  medium  having  parameters  Vp  =  4500 
m/Bec  and  y  =  1/ As  follows  from  Formula  (2.1)  and  (2.2),  the 
density  of  the  medlvim  evidently  does  not  enter  the  formula  for  the 
displacements  anywhere,  so  that  Vp  and  1  fully  describe  the  required 
mechanical  properties  of  a  medium  having  parameters  close  to  that  of 
rock  of  medium  hardness. 

Let  us  take  the  distance  30  for  R^,  where  R^  Is  the  radius  of 
the  charge;  the  depths  of  the  charges  are  20  R^  and  40  R^.  For  signal 
f(t),  used  In  the  calcvilatlon  and  having  a  duration  t*  »4,  some 
nominal  wavelength  may  be  taken  as  equal  to  X  =  RQt*  »  4  Rq  =  120  R^. 

This  also  determined  the  range  of  calculated  dimensionless  dis¬ 
tances  r  =  6o  to  240  R_. 

Let  us  examine  first  of  all  the  theoretical  oscillograms  of  dis¬ 
placements  at  a  distance  r  =  120  R,  for  h  -  20  R_  (Fig.  6)  .  The 
theoretical  oscillograms  of  r  =  120  R,,  h  =  20  R,  are  continuous  lines: 
1 — the  horizontal  component  of  motion,  and  2 — the  vertical.  In  the 
same  scale  along  the  abscissa  Fig.  6  shows  the  form  of  the  signal 
f(t)  (curve  3)  In  an  arbitrary  vertical  scale  (the  time  of  reading 
the  signal  Is  shifted  to  the  entry  time  of  the  wave  P^) .  In  com¬ 
parison  with  the  signal  f(t)  the  shape  of  the  curves  of  and 

w_  (t*)  Is  more  complex:  the  total  duration  of  the  signal  was  In- 
creased;  the  ratio  of  the  amplitudes  of  successive  peaks  and  troughs 
is  different  on  various  components;  and  Instead  of  a  simple  radial 
motion  In  the  wave  Pq  In  an  unbounded  medlvun  the  motion  of  a  siO'face 
particle  has  polarization  (Fig.  7) .  The  study  of  the  trajectory  of 
a  particle  Jointly  with  accurate  oscillograms  and  the  qualitative 
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pattern  of  §1,  Section  5  allows  It  to  be  concluded  that  various  waves 
take  part  In  the  complex  motion  of  a  surface  particle:  Pq  + 

+(Sj)  and  R.  If  we  turn  to  Pig.  1,  where  the  wave  fronts  Pq,  P^, 
eind  the  waves  Sj  and  R  are  shown  schematically,  we  can  Interpret  the 
complex  motion  of  the  surface  In  Pig.  7  as  follows.  The  beginning  of 
the  motion  of  a  svirface  point-above  and  ahead  of  the  epicenter — 
corresponds  mainly  to  motion  In  the  longitudinal  wave  P^  and  Includes 
reflection  effects:  Pj^,  S-^.  This  motion  Is  accomplished  elliptic  orbit 
with  clockwise  particle  rotation.  After  the  transitional  stage  with 
particle  motion  below  and  behind,  and  then  above  and  ahead,  particle 
motion  begins  In  a  counterclockwise  elliptic  orbit.  This  motion  Is 
Identified  with  the  wave  R,  and  the  transitional  stage  Is  Identified 
with  motion  In  transverse  waves  S  whose  envelope  Is  also  the  trans- 
verse  wave  S^. 


Pig.  7.  Trajectory  of  a  particle  of  the  surface  of  a  medium, 

according  to  data  In  Pig.  6. 

In  conclusion  to  the  qualitative  analysis.  It  should  be  noted 
that  separation  of  complex  motion  Into  Indlvldvial  waves  Is  conditional 
and  contains  the  following  concept.  Sepsu?ated  sections  of  theoretical 
oscillograms  and  trajectories  are  Identified  with  Individual  waves  on 
the  basis  of  kinematic  and  dynamic  properties  according  to  the  pre¬ 
eminent  manifestation  of  the  properties  of  the  Individual  waves  In 
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These  sections.  Generally  speaking  it  is  impossible  to  deny  the 
appearance  of  longitudinal  waves  in  section  Sj  and  of  trauisverse 
waves  Sj  in  section  R,  but  this  manifestation  of  additional  disturb¬ 
ances  is  weaker  than  the  fundamental  motions^  the  indices  of  which 
are  appropriate  to  the  section. 


Fig.  8a.  Theoretical  oscillograms  for  h  ■  20  R-*,  the  arrival 
time  of  wave  Pq  at  a  given  point  is  -^plotted  along 
the  abclssa  In^unlts  of  t*.  Symbols  of  components 
are  the  same  as  in  Pig,  6. 

Pig.  8b.  Theoretical  oscillograms  for  h  =  40  R^. 

Figs.  8a  and  8b  show  combined  theoretical  oscillograms  of  the 
motion  of  a  surface  for  these  same  parameters  of  the  mediiun^  for  the 
distances  and  depths  indicated  there. 
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Even  a  cursory  examination  of  the  oscillograms  In  Pig,  8a 
(h  =  20  Rj)  makes  It  possible  to  clearly  distinguish  two  fundamental 
groups  In  the  motion  of  points  on  the  surface  z  =  0,  which.  In  ac¬ 
cordance  with  the  qualitative  analysis,  are  denoted  by  the  Indices 
P  and  R  (Pq,  and  In  the  first  groiyp,  and  Sj  and  naturally  R 
In  the  second) .  These  groups  are  propagated  at  various  rates:  the 
first  at  v_  and  the  second  at  v„.  The  maximum  amplitudes  In  these 

P  A 

groups  decrease  with  an  Increase  In  r  according  to  various  laws: 

while  at  close  distances  6o  R^  the  first  group  predominates  In 

amplitude,  then,  starting  from  the  distance  ^  120  R_,  the  second  group 

Z> 

starts  to  predominate.  At  this  same  distance  motion  In  the  transi¬ 
tional  section  possesses  a  comparatively  Insignificant  amplitude, 
manifesting  Itself  mainly  In  the  phase  of  retrogressive  motion  In 
the  first  group  and  In  sections  of  growth  of  the  second. 


Pig.  9*  Theoretical  osclllo- 
greuns  for  r  =  120  R3,  h  *  20  R, 
emdir  »  0.4,  0.5  "^and  0.6.  ^ 
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On  the  basis  of  the  calculations.  In  spite  of  small  number  of 
points  with  respect  to  the  distances.  It  can  be  shown  that  at  close 
distances  the  maximum  displacements  In  the  first  group  are  decreased 
with  distance  not  slower  than  r"^,  while  at  the  end  of  the  calcvilated 
range  transition  is  projected  to  the  law  (I.15):  displacements  de¬ 
crease  more  rapidly  than  r”^,  but  more  slowly  than  r"^.  If  a  graph 
Is  constructed  of  the  damping  of  the  maximum  amplitides  of  displace¬ 
ments  or  velocities  according  to  1.15)  for  a  chosen  network  of  dis¬ 
tances  and  depth  h  =  20  R^,  It  is  possible  to  approximate  this  de¬ 
pendence  by  an  exponential  of  the  form  Ar“^,  where  n  =  1.6  to  1.7. 

This  law  is  the  same  for  both  displacements  and  particle  velocities. 

It  should  be  noted  that  It  becomes  applicable  earlier  (at  closer 
distances)  for  describing  particle  velocities.  In  the  calculated 
range  In  the  second  group  the  displacements  decrease  more  slowly  than 
r"^,  but  more  rapidly  than  r"^/^  and  more  closely  to  the  latter  law, 
which  Is  asymptotic  for  the  surface  wave  R.  This  comparative  charac¬ 
teristic  allows  change  of  the  carrier  of  maximum  amplltlde  to  be  In¬ 
dicated  In  general  motion  with  an  Increase  In  distance  from  the  epi¬ 
center.  In  this  the  first  group  may  be  traced  more  clearly  on  the 
horizontal  component  (first  maximum  and  first  minimum),  and  the  second 
wave  group  on  the  vertical  component,  at  the  end  of  recording. 

Fig.  8b  shows  theoretical  oscillograms  for  a  charge  depth  h  = 

=  40  R^.  A  comparison  of  the  oscillograms  In  Pigs.  8a  and  8b  In¬ 
dicates  that  the  group  which  Is  Identified  by  a  surface  wave  has  con¬ 
siderably  less  amplltlde  In  Fig,  8b,  although  some  decrease  In  ampli¬ 
tude  Is  noted  In  the  first  wave-group.  At  h  =  40  R^  change  of  the 
carrier  of  maximum  amplitude  does  not  take  place  In  the  entire  calcu¬ 
lated  range  of  distances.  From  a  comparison  of  the  damping  of  maximum 
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amplitudes  In  the  first  and  second  groups  It  can  be  proved  that  this 
change  at  h  =  4o  R_  also  takes  place,  but  at  greater  distances  than 
at  a  source  depth  h  =  20  Rj. 

In  conclusion  let  us  note  that  the  fundeunental  properties  of 
the  wave  pattern  do  not  In  essence  depend  upon  the  concrete  function 
f(t)  In  the  calculations.  The  only  Important  fact  here  Is  that  f(t) 
determines  a  certain  wavelength  X,  which  Is  the  characteristic  dimen¬ 
sion  In  the  problem,  and  this.  In  turn,  determines  the  range  of  dis¬ 
tances  r  »  X  under  study. 

3.  Let  us  examine  the  change  In  the  wave  pattern  occurring  with 
a  change  In  the  value  of  tt  .  The  parameter  v^  =  4500  nv'sec  Is  fixed. 
Change  In  7  In  this  case  denotes  a  change  In  the  propagation  rate 
of  transverse  and  svjrface  waves:  when  7  varies  from  0.6  to  0.4,  the 
value  Vg  varies  from  2700  m/sec  to  18OO  m/aec.  As  follows  from  the 
theoretical  oscillograms  (Plg.  9),  the  second  group  lags  behind  the 
first  more  as  7  decreases.  In  this  It  Is  possible  to  outline  (7  = 

*  0.5)  and  clearly  separate  (7  =  0.4)  the  Intermediate  maximum  In 
motion  In  the  transitional  section,  which  corresponds  to  the  second 
maxlmvun  of  the  signal  f(t)  .  All  fundamental  properties  of  the  wave 
pattern  noted  In  Section  2  are  preserved  when  7  varies.  The  maxlmvun 
amplitudes  of  displacements  when  7  varies  from  0.6  to  0.4  differ 
negligibly  from  the  general  law  of  dependence  Indicated  by  the  scale 
factor  k  =  U()R(y^Vg.  This  follows  from  the  graphs  In  Pig.  9,  where  all 
curves  are  constructed  In  a  single  vertical  scale  k  =  uqRo/iSoo. 

In  conclusion  let  us  point  out  that  this  qualitative  analysis  of 
the  wave  pattern  of  the  motion  of  a  free  surface  of  an  elastic  medlvim 
should  be  supplemented  by  an  analysis  of  the  changes  In  spectrvun  and 

energy  In  elastic  motion  of  the  boundary.  These  problems  are  separ¬ 
ate  due  to  their  complexity  and  size  and  will  be  the  subject  of 
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special  discussions. 

The  Author  thanks  K.  Ye.  Oubkln,  Y.  N.  Rodionov  and  A.  N.  Romashov 
for  valuable  discussion  of  the  results  of  work  In  the  preparation  of 
this  article. 
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DiMOLITION  OP  A  SOLID  BY  AN  EXPLOSION 
V.  V.  Adushkln  and  A.  P.  Sukhotln 


In  this  work  we  will  examine  the  demolition  of 
a  Plexiglas  block  by  the  explosion  of  a  spherical 
charge  Inside  It.  The  motion  occurring  during  the 
explosion  has  a  spherical  symmetry  which  facilitates 
observation  and  Interpretation  of  the  experimental 
results.  We  will  note  that  to  study  the  qualitative 
picture  of  demolition  the  selection  of  the  solid  be¬ 
ing  demolished  by  the  explosion  Is  not  Important. 

The  experimental  data  cited  In  this  work  made  It 
possible  to  elucidate  certain  aspects  of  the  de¬ 
molition  of  a  solid  by  an  explosion.  A  general 
plctvire  of  destruction  Is  obtained  and  the  time 
characteristics  of  the  process  are  evaluated.  It 
Is  shown  that  during  demolition  (both  In  the  Inside 
zone  and  with  the  shatter)  It  Is  necessary  to  dis¬ 
tinguish  two  aspects:  the  conception  of  cracks  and 
their  development.  It  turned  out  that  the  space- 
time  characteristics  of  the  demolition  process  sub¬ 
stantially  depend  on  the  behavior  of  the  material  In 
the  Immediate  vicinity  of  the  explosion  point.  The 
presence  of  a  plastic  zone  near  the  discharge  Is 
equivalent  to  an  Increase  In  the  dimensions  of  the 
source  radiating  an  elastic  wave. 
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Description  of  experiments.  Plexlglas(poly( methyl) methacrylate) was 

selected  as  the  solid  to  be  demolished  by  explosion.  We  will  cite  the 
physical  characteristics  of  Plexiglas  taken  from  the  handbook  [l]: 

Density  .  I.l6 — 1.20  g/cm^ 

Sound  velocity .  2800 — 5000  xa/sec 

Heat  capacity . 0.55 — 0.40  cal/g  •  deg 

Thermal  conductivity . (l — 10)  *10"^  cal/cm*  sec ‘deg 

Tensile  strength . -^1000  kg/cm^ 

The  blocks  were  cemented  together  from  Plexiglas  sheets  50-50  mm 
thick  with  V-51-P9  glue  or  with  a  5  %  solution  of  Plexiglas  In  dlchloro- 
ethane.  Several  experiments  were  made  with  monolithic  blocks.  We  will 
cite  the  dlmenslcxis  of  the  blocks  which  were  tested: 
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In  the  center  of  the  block  we  placed  a  charge  of  pressed  PETTN 
weighing  0.8  g  with  a  radius  of  r^  =  5  nun,  or  a  cast  charge  of  TNT- 
RDX  50/50  weighing  24  g  with  a  radius  of  r^  =  15  mm.  In  the  experi¬ 
ments  with  monolithic  blocks,  the  charge  was  placed  In  a  hole  vrtilch 
was  closed  by  a  threaded  Plexiglas  stopper. 

The  development  of  the  explosion  and  demolition  of  the  block  was 
photographed  In  the  passage  of  light  by  means  of  a  high-speed  camera 
(SFR),  which  In  addition  to  continuous  time  scanning  permitted  us  to 
obtain  a  series  of  frames  with  a  frequency  of  up  to  2  million  frames 
per  second.  Plgxore  1  shows  the  diagram  for  carrying  out  the  experiments. 
Block  2  was  set  so  that  the  siu’face  of  the  splice  was  parallel  to  the 
light  rays  coming  from  the  pulse  laiup  (lPK-120)  through  condenser  lens 
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3.  Oetpnatlon  of  charge  1  of  the  explosive  was  done  by  a  high-voltage 
pulse  obtained  by  discharge  of  condenser  KBO-P  (C  =  6  Mf  )  charged 
to  a  voltage  of  4  kv.  Breakdown  of  the  discharge  gap  In  the  detona¬ 
tion  circuit  was  Initiated  by  an  electrical  pulse  supplied  from  con¬ 
trol  panel  4  (SFR)  .  The  Instant  of  Its  delivery  was  coordinated  with 
the  position  of  the  scanning  mirror.  This  pulse  switched  on  the  boost¬ 
er  lamp. 

Results  of  the  experiments.  Figure  2  shows  typical  film  frames 
of  the  demolition  of  a  100  x  100  mm  block.  Figures  3  and  4  show  the 
continuous  time-scanning  of  this  process  in  blocks  of  different  size 
(W  denotes  the  distance  from  the  center  of  the  charge  to  the  free 
surface) .  The  wavefront  and  demolition  zone  are  distinctly  seen  In 
the  photographs.  The  general  picture  of  the  development  of  the  basic 
phenomena  during  explosion  of  a  PETN  charge  weighing  0.8  g  In  a  Plexi¬ 
glas  block  Is  shown  In  r-t  coordinates  In  Pig.  5*  Line  I  represents 
the  track  of  the  compression  wavefront  which  Is  reflected  from  the  free 
surface  as  a  tension  wave  (line  2) .  The  contour  of  region  3  shows 
the  expansion  of  the  charge  chamber  In  time.  Its  development  ceases 
In  approximately  6  m  sec  when  the  radius  reaches  the  value  Pr^.  Start¬ 
ing  from  distance  3rQ,  about  6  ft  sec  after  the  passage  of  the  compres¬ 
sion  wavefront  there  Is  a  new  region  In  which  the  material  Is  demolish¬ 
ed,  region  4.  Line  4  on  the  graph  shows  the  development  of  Its  dimen¬ 
sions  In  time.  We  see  from  the  graph  In  Pig.  5  that  In  a  block  100  mm 
thick  (20r^)  the  effect  of  the  free  sxirface  (line  2)  on  the  develop¬ 
ment  of  region  4  Is  absent  at  least  during  24  nsec. 

Compression  wave.  The  graph  In  Pig.  5  and  the  photographs  show 
that  the  velocity  of  the  ccngjresslon  wavefront  varied  In  the  range 
of  distances  to  (5-6) r^.  The  sample  graph  of  the  dependence  of  the 
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velocity  of  the  front  on  the  distance^  which  was  obtained  by  a  graphic 
differentiation  of  line  l,  is  shown  In  Fig.  6.  The  dashed  line  in 
the  graph  marks  the  value  of  the  velocity  of  the  longitudinal  acoustic 
waves  In  the  Plexiglas.  We  see  from  the  graph  that  In  the  range  from 
Ir^  to  (5-6)rQ,  the  velocity  of  the  compression  wavefront  varies  from 
5000-6000  nv^sec  to  2900  m/aec.  The  velocity  of  the  front  of  the  re¬ 
flected  wave  (line  2  in  Pig.  5)  Is  somewhat  slower  than  the  velocity 
of  the  longitudinal  acoustic  waves  and  Is  less  than  2600-2700  m/aec. 

In  some  cases  we  see  parallel  to  line  1  still  another  line  (Plg. 

7)  whose  origin  Is  unexplained. 

The  magnitudes  of  the  mass  velocities  in  the  compression  wave  can 

be  estimated  by  the  Initial  velocities  of  the  movement  of  the  free 

svu*face  Vq>  assvunlng  that  the  velocity  doubles  with  reflection.  The 

Initial  velocities  of  the  motion  of  the  free  surface  at  the  instant  of 

reflection  of  the  compression  wave  were  measured  within  about  10%. 

We  will  give  the  values  of  the  mass  velocities  U  =  Vq/2  meters  per 

second,  as  well  as  estimates  of  the  radial  stresses  In  the  wavefront 

calculated  by  the  formula  4 _  =  cU: 

•  rr  o 

r/r«=|  5  6.8  7.',  JO 

=  100  75  55  30  .«/sec 

=J3500  2600  1900  1050  Icg/CM* 

If  we  assume  that  at  the  Interface  with  the  charge  the  pressure 
In  the  wave  eqxials  10^  kg/cra^,  we  can  write  approximately  the  law  of 
attenuation  of  the  maximum  pressure  In  the  conpresslon  wave  with 
distance  as 
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])ffmQllfelon  of  block.  An  examination  of  the  pieces  of  the  smashed 
block  permitted  us  to  refine  certain  details  of  the  picture  seen  on 
the  film  frames.  In  all  experiments  with  Plexiglas  blocks,  regardless 
of  their  dimensions,  we  observed  the  following  typical  picture  of  de¬ 
molition.  The  Plexiglas  In  a  range  of  distances  from  Ir^  to  Jr^  has 
no  apparent  destruction,  remaining  transparent.  In  the  vicinity  of 
the  charge  there  takes  place  a  noticeable  heating  of  the  material 
(At*  lO^C)  .  As  the  evaluations  showed  It  Is  Impossible  to  explain 
the  rise  In  temperature  as  the  result  of  heat  conduction.  It  Is  neces¬ 
sary,  apparently,  to  assume  that  the  material  Is  heated  due  to  shock 
compression  and  plastic  deformation  behind  the  wavefront.  At  distances 
of  we  observed  a  spherical  zone  consisting  of  a  large  amoxint 

of  fine,  radially  oriented  cracks.  The  number  of  these  cracks  per 
unit  surface  decreases  with  distance.  The  orientation  of  the  cracks 
attests  to  the  fact  that  they  developed  as  the  resiilt  of  the  action  of 
azimuthal  tensile  stresses. 

It  Is  Interesting  to  note  that  the  destruction  occurring  during 
reflection  of  a  compression  wave  from  the  free  surface  ( shatter)  has 
exactly  the  same  character.  The  resemblance  of  the  form  of  the  cracks 
In  the  shattered  and  central  zones  confirm  that  tensile  stresses  cause 
radial  cracks  In  the  range  (5-^)1*^.  We  already  noted  (graph  In  Pig. 

5,  range  4)  that  the  radial  cracks  In  range  (5-4)rQ  develop  about 
6  Msec  after  the  passage  of  the  compression  wavefront.  This  means 
that  azimuthal  tensile  stresses  are  absent  In  the  compression  wave 
front.  They  occur  later  as  a  result  of  the  radial  displacement  of 
the  material  In  the  spherically  spreading  compression  wave.  Figure 
8  shows  the  dependence  of  the  velocity  of  the  front  of  radial  cracks 
on  the  distance  which  was  obtained  by  a  graphic  differentiation  of 


curve  4  on  the  griiph  In  Pig.  5.  It  follows  from  the  graph  In  Pig.  8 
that  the  front  of  the  radial  cracks  after  Its  development,  at  a  dis¬ 
tance  of  5r^  has  a  velocity  of  about  2000  ra/sec,  which  decreases  to 
1000  va/Bec  at  a  distance  of  4r^.  The  magnitude  of  the  velocity  In 
this  region  Is  probably  determined  by  the  propagation  of  a  certain 
phase  of  aslmuthal  tensile  stresses  and  not  by  the  velocity  proper  of 
the  development  of  cracks.  This  Is  particularly  Indicated  by  the  fact 
that  the  fractvire  zone  consists  of  fine,  separate  cracks  not  connected 
together  In  the  direction  of  propagation. 

When  r  >  4r^  we  observe  a  change  In  the  character  of  destruction; 
only  Individual  cracks  develop  farther,  forming  a  wavy  type  of  fracture. 
Their  development  takes  place  with  a  velocity  decreasing  from  1000 
n^sec  to  500  nv^sec  at  a  distance  of  (5-6)rQ.  With  a  farther  propaga¬ 
tion,  the  velocity  of  the  cracks  does  not  change.  Apparently  In  the 
range  exceeding  4rQ  the  tensile  stresses  In  the  wave  attenuate  so  much 
that  they  become  Insvifflclent  for  generating  new  cracks. 

The  time  role  of  the  effect  of  the  explosion  products  compressed 
In  the  cavity  can  be  qvialitatlvely  exposed  from  a  comparison  with  the 
experiment  where  the  explosion  of  a  charge  occvirred  In  a  hole  on  the 
free  surface  of  a  plate  of  Plexiglas  50  mm  thick  ( lOr^) .  The  result 
of  the  experiment  Is  plotted  on  the  graph  In  Pig.  5^  where  curve  ^2 
shows  the  development  In  time  of  the  zone  of  radial  cracks,  region  4. 
The  picture  of  the  development  of  an  explosion  in  this  case  was  ob¬ 
tained  Just  as  In  the  case  of  an  explosion  Inside  the  block.  A  sub¬ 
stantial  difference  Is  that  the  single  radial  cracks  cease  growing 
already  at  a  distance  of  (4-5) r^.  This  can  be  explained  by  the  fact 
that  In  the  case  of  an  explosion  on  the  surface,  the  explosion  pro¬ 
ducts  have  the  opportunity  to  spread  out  freely  In  the  atmosphere  as 
a  result  of  which  the  stresses  In  the  material  are  reduced. 


Pig.  2.  Photographs  of  the  development  of  the  explosion  of 
a  PETN  charge,  C  «  0.8  g.  In  glued  block  of  Plexiglas.  The 
time  passing  from  the  Instant  of  explosion  Is  shown  under 
the  frames  In  microseconds. 
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Pig.  5*  Explosion  of  TNT-RDX 
50/50,  C  =  24  g.  In  Plexiglas 
block  (2W  =  lOr^ 


Fig.  5.  r-t  diagram  of  the 
development  6f  the  explosion 
In  a  block  with  thickness  2W= 

=  SOTq. 


Fig.  7.  Explosion  of  P; 
in  Plexiglas  block. 


Pig.  4.  Explosion  of  PETN 
charge,  C  =  0.8  g.  In  Plexi¬ 
glas  block  (2W  =  20rQ)  . 


Pig.  6.  Velocity  of  compression 
wavefront  vs.  distance. 


Fig.  8.  Propagation  velocity  of  the  zone 
of  radial  cracks  vs.  distance. 


It  is  interesting  to  note  the  effect  of  gluing  on  the  picture  of 
demolition.  The  film  frames  in  Pig.  2  refer  to  a  block  glued  by  di- 
chloroethane .  It  is  apparent  that  gluing  disrupts  the  properties  of 
the  Plexiglas,  forming  relieving  swfaces.  The  experiments  with  blocks 
glued  by  V-51-P9  glue  show  that  the  semlspherlcal  regions  of  shatter¬ 
ing  are  close  together  and,  finally,  the  results  of  the  experiment  in 
the  monolithic  block  attest  to  the  fact  that  if  the  properties  of  the 
material  are  the  seune  in  all  directions,  the  zone  of  destruction  is  a 
spherical  surface.  Figure  9  shows  the  film  frames  of  this  experiment 
with  a  120  X  120  mm  block. 

Expansion  of  the  discharge  chamber.  The  development  of  a  spherical 
dark  spot,  a  cavity,  (Fig.  2)  is  observed  in  the  center  of  the  block 
at  the  explosion  site  after  the  flash.  Region  5  represents  the  cavity 
on  the  graph  In  Pig.  5.  The  dependence  of  the  radius  of  the  cavity 
rj^  on  time  Is  shown  In  Pig.  10.  The  corresponding  dependences  for  a 
gas  bubble  In  water  and  a  cavity  in  sand  are  shown  on  the  same  graph 
for  comparison. 
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It  Is  Interesting  that  for  the  indicated  media  approximately  the 
same  law  of  motion  was  obtained  in  the  initial  expansion  section 


In  the  case  of  Plexiglas  the  maximum  redlus  of  the  cavity  based 
on  the  photographs  reaches  a  magnitude  of  Pr^. 

The  measvirements  of  the  radius  of  the  cavity  where  the  explosive 
had  been,  based  on  the  fragments  of  the  blocks,  show  that  the  cavity 
radius  Increases  on  the  average  by  only  20% (an  expansion  of  the  cavity 
to  1.2rQ) .  A  special  test  showed  that  optical  jffects  here  are  ex¬ 
cluded.  At  first  glance  such  a  contradictory  behavior  of  the  cavity 
can,  app2irently,  be  explained  by  the  experimentally  known  capacity 
of  Plexiglas  for  large  reversible  deformations  during  heating.  Under 
the  effect  of  the  compression  waves  and  pressure  of  the  explosion 
products,  the  heated  Plexiglas  is  expanded  to  the  size  shown  In  the 
photographs.  The  destruction  of  the  block  releases  the  compressed 
explosion  products.  This  results  In  the  elastic  forces  of  the  material 
compressing  the  cavity  to  a  size  of  1.2r^. 

Shattering.  Destruction  In  the  vicinity  of  the  free  surface  of 
the  block  (shatter)  Is  caused  by  tensile  stress  which  arises  during 
the  reflection  of  a  wave  from  the' free  surface.  Destruction  of  this 
type,  which  Is  observed  in  the  present  experiments,  occurred  practi¬ 
cally  at  the  front  of  the  reflection  wave.  An  examination  of  the 
fragments  after  the  experiment  showed  that  the  shatter  crack  In  the 
Plexiglas  block  Is  not  a  single  continuous  crack  but  Is  a  certain  zone 
of  destroyed  material  ^-5  mm  thick  consisting  of  numerous  small  cracks 
In  the  form  of  lentlcles  or  scales  (0.5-3. 0  mm  In  size) .  The  shatter 
cracks  are  oriented  parallel  to  the  front  of  the  reflection  wave  and 
only  as  a  result  of  their  development  form  an  Interface.  When  the 
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ireseyve  of  kinetic  energy  Is  siifflclent  In  the  piece  Intended  for 
separation,  the  Individual  cracks  of  the  boundary  of  the  shatter  zone 
grow,  reach  the  free  surface,  and  the  small  shatter  plate  files  off. 

The  noted  type  of  shatter  destruction  and  the  fact  that  the  shatter 
cracks  which  appeared  are  developed  In  time  compels  us  to  re-examine 
the  frequently  used  assumption  of  the  Instantaneous  destruction  of 
material.  It  la  necessary,  apparently,  to  estimate  carefully  the 
time  of  destruction  In  each  specific  case.  In  the  assumption  of  In- 
staintaneous  destruction  the  head  portion  of  the  wave  pvilse  Is  coinplete- 
ly  cut  off  by  a  crack.  A  certain  delay  In  the  formation  of  a  shatter 
crack  caused  by  Its  development  can  result  In  the  head  part  of  the 
wave  cxirve  In  the  stress-time  coordinates  changing  Its  form  In  the 
region  where  the  shatter  cracks  are  developing.  This  In  turn  Inevitably 
entails  a  change  In  the  character  of  shatter  destruction  as  a  whole. 
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Fig.  10.  Dependence  of  the 
radius  of  the  cavity  on  time: 
1)  Plexiglas,  A  =  270;  2) 
water,  A  =  290;  5)  sand,  A  = 
=  250. 


Pig.  11.  Distribution  of  the 
radial  components  of  the  stress¬ 
es  behind  the  wavefront. 


)llcablllty  of  an  Elastic  Model  of  a  Medium  When  Analyzlr 


an  Explosion  In  a  Solid.  Let  us  consider  the  simplest  problem.  In 
an  elastic  unbounded  space  with  a  Poisson  ratio  •  «  l/5  In  a  cavity 
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with  radius  r-|^  at  the  Instant  of  time  t  =  0  there  Instantaneously 
arises  a  pressure  Pq  which  then  remains  constant.  Figures  11  and  12 
show  the  results  of  the  calculation  In  the  form  of  a  dependence 

<3»r  /  />o  =  /i  ('•  /  /'i),  3i,o  /  Po  =  /•:  ('•  /  '-i) 

at  different  Instants  of  time,  where  Is  the  stress  In  a  radial 
direction,  and  Is  the  stress  In  the  azimuthal  direction. 

Let  us  compare  the  length  of  the  positive  phase  of  the  compression 
wave  ,  estimated  by  the  thickness  of  the  shattered  plate  8,  with 
the  calculated  value  (Plg.  ll)  .  In  the  experiment  with  the  block 
with  the  size  2W  =  20r^  (rQ  =  5  nun)  the  thickness  of  the  shatters 
plate  was  10  mm.  In  this  case  the  wave  length  X^  -  2  8  =  20  mm, 
which  corresponds  to  ^r^. 

It  follows  from  the  solution  of  the  problem  that  for  this  case 
r  =  10  r^,  the  value  X^  =  I.^^q  =6.5  mm,  l.e..  In  the  experiment  X^ 

Is  considerably  greater  than  should  have  been  In  an  elastic  medium. 

An  examination  of  the  graph  In  Pig.  12  shows  that  large  tensile 

stresses  <r  are  observed  only  In  the  neighboring  zone  at  distances 
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of  (1.5-2) r^.  Beyond  the  limits  of  this  zone  sharply  drops.  The 
experiment  shows,  however,  that  the  central  zone  of  demolition  has  a 
size  up  to  (5«5-^)rQ.  Thus  a  comparison  of  the  experimental  data  and 
of  the  results  of  the  calculation  with  respect  to  the  elastic  model 
attest  to  the  fact  that  this  model  cannot  be  used  for  describing  the 
behavior  of  a  material  In  the  entire  range  of  distances. 

The  difference  between  the  experiment  and  the  calculation,  as  well 
as  the  character  of  the  observable  destruction  of  the  material  proclaim 
that  In  the  neighboring  zone  the  material  behaves  explicitly  Inelastic- 
ally.  Thus,  In  the  range  of  about  we  observe  a  substantial  heat¬ 

ing  of  the  material  which  Indicates  the  presence  In  this  zone  of  In- 
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elastic  deformations.  As  distance  from  the  charge  Increases  the  stres¬ 
ses  In  the  wave  drop  and  we  can  expect  that  with  a  certain  distance 
r  =  r^  the  elastic  model  of  the  medlvun  will  apparently  be  Justified 
to  some  extent.  We  will  try  to  coordinate  the  qualitative  picture, 
having  selected  as  an  appropriate  sample  a  radius  of  the  source  of 
elastic  oscillation  r^^.  If  we  select  as  r^^  the  radius  of  the  plastic 
region  2.5rQ,  then  we  will  obtain  a  satisfactory  agreement  of  the 
value  of  the  shatter  plate  and  the  calculated  wavelength  and  also  Of 
the  radius  of  the  formation  of  nuclei  of  radial  cracks  and  of  the 
zone  of  strong  azimuthal  stresses.  We  can  assume  that  the  plastic 
flow  of  material,  which  Is  observable  dvirlng  explosion,  leads  to  an 
Increase  in  the  effective  radius  of  the  source  In  an  elastic  model. 


Pig.  12.  Distribution  of  azimuthal  components  of 
stress  behind  the  wave  front. 
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